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Abstract 



The multipole amplitudes for the N - A(1232) electromagnetic transition are 

computed in the framework of the linear a model and the chiral chromodielec- 

•^ ' trie model for small and moderate photon virtualities. The models include 

Qv ' quark and meson degrees of freedom and the nucleon and the delta are clusters 

^P . of three valence hedgehog quarks surrounded by meson clouds described by 

coherent states. Angular momentum and isospin projections are performed to 

endow model states representing the nucleon and the delta with proper quan- 

Q_i' tum numbers. Recoil corrections involved in the process ^^N -^ A are taken 

Q^' into account by performing linear momentum projection of the initial and fi- 

D . nal baryon states. The ratios E2/M1 and C2/M1 are in good agreement with 

the data in the two models, but the magnetic amplitude is better reproduced 

in the Linear Sigma Model. The ratios show little dependence with the model 

/\ ' parameters. Both in the Linear Sigma Model and in the Chromodielectric 

C^ ' Model the charged pions are responsible for the non-vanishing quadrupole- 

electric and -coulomb amplitudes. The recoil corrections enhance the results 

obtained for the amplitudes without linear momentum projection, improving 

the comparison with experimental data. The dependence of the theoretical 

amplitudes with the choice of the reference frame is also studied. 
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I. INTRODUCTION 

Electromagnetic processes have always played a central role in studies of the structure 
of nuclei, nucleons and their excitations. Recently, the interest in the electroproduction 
of the A(1232) and other nucleon resonances has increased, fueled by the large number of 
experiments planned and already running in several centers (Mainz, Bonn, MIT, TJNAF, 
etc) , where very clean electromagnetic probes are now available. In this paper we report on 
a theoretical calculation of the multipole amplitudes of the N-A electromagnetic transitions. 

The process 'j^N — > A has been considered in the framework of several models of baryon 
structure. From the point of view of a pure quark model, the A state results from a spin flip 
of one quark in the nucleon. This corresponds to a magnetic-dipole transition and vanishing 
quadrupole transitions. The experimental observation of quadrupole transitions, although 
small in comparison with the magnetic dipole, caused a discussion about the structure of 
the nucleon and the delta. In models only with quarks, non-vanishing quadrupole electric 
and scalar nucleon-delta transition amplitudes result from d-state admixtures to the quarks' 
lowest s-state, otherwise those amplitudes would be identically zero |jl]-§|. The quadrupole 
transitions resulting from such charge deformations are generally small. However, other 
explanations can be found, in particular the contribution of pions included in the baryon 
model states 1^-0 , as in the type of effective theories considered in the present work. 

Our calculations are carried on in the framework of two well known quark-meson mod- 
els of the baryon structure, namely, the Linear Sigma Model (LSM) |^ and the Chiral 
Chromodielectric Model (CDM) ||^, which have been used to describe the structure of the 
nucleon ^-Q. In these models, a baryon — such as the nucleon — is a soliton with three 
bare valence quarks, all in the same orbital-spin-isospin state, interacting with chiral a and 
7? meson fields. In the CDM, there is an additional interaction with a scalar-isoscalar chiral 
singlet meson field — the chromodielectric field, x- Except for a small explicit chiral sym- 
metry breaking term, both models are SU(2)xSU(2) chiral invariant, a symmetry which is 
spontaneously broken to SU(2), the pions being the Goldstone bosons. 

Although the two models use essentially the same ingredients, they provide quite differ- 
ent pictures of the nucleon (and of the delta). In the LSM the stability of the cluster of 
three quarks interacting with the mesons depends on the quark-meson interaction strength. 
A soliton is formed when the coupling constant for that interaction is sufficiently large, and 
it turns out that a strong meson cloud (particularly a pion cloud) is required for stabilizing 
the system. In the LSM, the chiral mesons bind the three quarks. In the CDM, besides 
the interaction between the quarks and the chiral mesons, there is an interaction between 
the chromodielectric field and the quarks. As a result of this interaction the quarks acquire 
a position-dependent dynamical mass which is an increasing function of the distance. The 
quarks are thus prevented to move too far away from the origin and such mechanism effec- 
tively generates quark confinement. The role of the chiral mesons, although conceptually 
important to implement chiral symmetry and its dynamical breaking, gets much suppressed 
and the resulting picture of the baryon is a soliton with three confined valence quarks sur- 
rounded by a weak cloud of chiral mesons (particularly pions) . By considering two models 
providing such extreme pictures of the nucleon, namely with quite different meson clouds, 
we are able to address interesting questions e.g. how electromagnetic transition amplitudes 
are sensitive to quantities such as the "number of pions" presented in the baryon states. 



In the framework of the chiral models considered in this paper, the nucleon and the 
delta are made out of neutral {a, x and 7r°) and charged (quarks and vr^) particles, the 
latter coupling directly to the electromagnetic probe (a virtual or a real photon). The three 
quarks are assumed in the same orbital-spin-isospin state occupying the lowest s-state (no 
(i-state admixture), represented by |g). Moreover the hedgehog state is assumed for the 
spin-isospin wavefunction of the quarks. The three valence quarks are therefore described 
by the fully symmetric state, |g^) (anti-symmetrization applies in color space — the state 
is a Slater determinant in color space). A quantum mechanical description for the mesons 
is considered by means of coherent states representing pion, sigma and chi clouds, namely 
|n), |S) and \x)- The starting point to describe a baryon in the framework of the LSM and 
CDM is, therefore, the Fock state \ip) = \q) |n) |S) \x)- Such state should then be projected 
onto angular momentum and isospin eigenstates in order to get states with the nucleon and 
delta quantum numbers ||lOHl2|. 

The calculation reported in this work is a natural extension of [Q. We have refined the 
approximations, namely by taking into account the state of motion of the initial and final 
baryon states involved in the nucleon-delta transition. To this end, a linear momentum 



projection [^ of the initial and final baryon states is applied, following the method used 
in []T1| for the calculation of the nucleon electromagnetic form factors. Besides this concep- 
tual improvement, we also present in more detail the formalism and address the issue of how 
the choice of the reference frame affects the theoretical transition amplitudes. 

Other calculations of the electromagnetic A^-A transition amplitudes carried out in sev- 
eral effective models of the nucleon have been reported in constituent quark models |l|-§I, 



with two-body exchange currents |T^, in the Skyrme model [p!6| , p!7| , in the cloudy bag model 



0J^, in chiral quark solitons of Nambu-Jona-Lasinio type with polarized Dirac sea |jT8|,|T9 
etc. 

This paper is organized as follows. In Section |I| we give a short account of the models 
and sketch the approximations used to construct model states representing the nucleon and 
the delta. In Section |T| we develop the formalism for the application of the models to 
the electroproduction of the A(1232) with a special emphasis on the implementation of the 
linear momentum projection. Finally, the results are presented in Section |IV| together with 
their discussion. A summary of the main conclusions of this work is presented in Section 0. 
The more technical aspects are given in the appendices. 

II. MODELS AND MODEL REPRESENTATION OF BARYONS 

The Lagrangian densities of the LSM and the CDM can be written, in a compact form, 
as 

where 

/:g = v^l^d^^ (2) 

li,,. = \{daf ^]^{d^f -V{a,ii) (3) 



are the pure quark and chiral meson terms, 



A-^,^,X = — V^ (CT + if ■ 7f75) i) (4) 

A. 



is the quark-meson interaction term and 
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^x = o(5x)^ - -Mlx\ (5) 



absent in the LSM, contains the kinetic and potential terms for the chromodielectric field. In 
these expressions %Ij{x) represents the quark field operator, ti^x) and a{x) the chiral pion and 
sigma meson fields, respectively (the arrow denotes isovector), and x(x) the chromodielectric 
field. The parameter p in the denominator of the interaction Lagrangian £g_o-,7r,x is in the 
LSM (no X fidd in this model) and 1 in the CDM. 

The meaning of the other terms appearing in eqs. (0)-(^) is the following. In Eq. (^, 
t/((T, vr) is the Mexican hat potential, g in Eq. (^) is the coupling constant which is dimen- 
sionless in the LSM and has dimensions of energy in the CDM (with p = 1). In Eq. (^ the 
second term on the r.h.s. is just the mass term for the x field, M^ being its mass. Other 
versions of the CDM consider a potential which includes, besides the mass term, up to quar- 
tic terms in the x field, as well as other powers of p in the interaction term (^. By just 
taking the mass term in the potential for x and p = 1 m the interaction, quark confinement 
is imposed in the smoothest way, which is the most appropriate choice for the quark matter 
sector of the CDM [^ . The Mexican-hat potential is given by: 



U = ^(a^ + T^^-v^f + ca + d. (6) 

The SU{2) X SU{2) chiral symmetry of C is explicitly broken by the small term ca. The 
parameters A, z/ and c are related to the sigma and pion masses, m„ and m^, and to the 
pion decay constant, /t^: 

ml -ml 2 .2 "^^ r 2 /7^ 

= — ^ — ' ^ ^ ^''~~' c = -Um^. (7) 

In (^, (i is a constant which guarantees that min U = 0. The Mexican hat potential induces 
spontaneous chiral symmetry breaking. The vacuum expectation values of the chiral fields 
are zero for the pion and — /j^ for the sigma: 

(0|#|0)=0 (8) 

{0\cr\0)=-U (9) 

(we use the hat symbol """ whenever we want to stress the operatorial character of the 
fields). It is convenient to define a new sigma field, which we still denote by a, as the 
fiuctuation around the vacuum value, i.e., we perform the replacement a —^ — /tt + o". Hence, 
the vacuum expectation value of the 'new' sigma field is zero, according to @. 

Altogether, the parameters of the models defined by (|T]) are the pion and sigma masses 
(fixed at m^ = 0.139 GeV and m.^- = 1.2 GeV), the pion decay constant (/^ = 0.093 GeV), 
and g in the LSM and g and M^ in the the CDM. In the simpler version of the CDM 



G 
G 



gM^. 
0.2 GeV, 



considered in this work, it turns out that the results are sensitive only to the combination 
The physical region of the coupling constant is (? ~ 5 in the LSM |pi]|JT^ and 
in the CDM |n,2^. We remark that the physical range of the coupling 
constant in the CDM is much narrower than in the LSM. For this reason, later on, in 
Section IV, when we show the dependence of the results on the coupling constant (all other 
parameters being fixed to the quoted values) we shall only consider the LSM. The only free 
parameters are g in the LSM and G in the CDM and these are fixed in order to reproduce 
well the bulk of the nucleon properties. 

It is known that, in these quark- meson models, the delta- nucleon mass splitting is small. 
We may remedy this by adding to the Hamiltonians of the models explicit bare baryon mass 
terms with different masses for the bare nucleon and the bare delta @,0. Then there is 
one more parameter — the bare nucleon-bare delta mass difference — which can be fitted 
to reproduce the physical nucleon-delta mass splitting. The inclusion of such a term has 
little effect (specially in the LSM) in the wave-functions of both quarks and mesons. Such a 
bare nucleon-delta mass splitting accounts for the residual chromomagnetic interaction and 
for the 't Hooft interaction which is attractive for the nucleon and absent in the delta. 

Solutions of the LSM and CDM representing the physical baryons can be obtained using 

]JT|j2gj2|]. For the sake 



a variational approach based on the projected hedgehog ansatz 

of completeness we sketch the formalism here. We consider three valence quarks with spin 

and isospin state in the so-called hedgehog configuration: 
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\hh) = -^{\ui) 



MT)). 



(10) 



All quarks occupy the same lowest positive energy s-state of the model effective potential, 
given by the spinor 

1 / u{r) 
icr ■ rvir) 



Qh) 



\hh) 



(11) 



In our approach, the pions, sigmas and chis (in the CDM) are described by coherent 
states: |n) for the pions, |S) for the sigmas and \x) for the chis. The expectation values of 
the field operators in these coherent states are the mean meson fields. The hedgehog ansatz 
for the mesons reads: 



(S|a(r)|S) =(j{r) 

(n|#(r)|n) = -0(r) 
r 

(xlx(^)lx) =x{r) 



(12) 
(13) 
(14) 



(we remember that a is now the fluctuating part of the original sigma field around the 
vacuum expectation value, —f-n) ■ Actually the spherical symmetry of sigmas and chis and the 
"hedgehog-like" character of the pion with the peculiar isospin-coordinate space correlation, 
result from the quark spin-isospin hedgehog configuration (|10|) and from the requirement of 
minimum mean field energy solutions ||2l|-^ 
The pion coherent state is 

\Il)=NMe^^\Y,fSk 



luj^{k) 



Uk)a\{k) |0) 



(15) 



where aj(fc) creates a free pion with momentum k and (Cartesian) isospin index i, N.^ is a 



normahzation factor, cu^ = Jk"^ + m^ and ^i{k) is the pion amphtude. Similarly, the sigma 
coherent state is given by 



|S) = N^ [r,] exp j d'ksj'^ r^{k)h\k) |0) , (16) 

where 6^(fc) is the sigma creation operator, r]{k) is the coherent state amplitude function for 
the cr-field. A similar expression holds for the chi field and we denote the amplitude of the 
corresponding coherent state by K,{k). 

The coherent states are particularly easy to deal with because they are eigenstates of the 
annihilation operators, e.g. 



a,(fc)|n) = ^^6(fe)|n). (17) 

Similar expressions hold, involving the annihilation operator of sigmas and chis. 

The coherent state amplitudes are the Fourier transforms of the meson functions in coor- 
dinate space introduced in (p!^)-(|T^), and exhibit the following hedgehog shape in momentum 
space: 

Uk) = -^^m (18) 

r]{k) = ri{k) (19) 

K{k) = K{k). (20) 

Altogether the hedgehog baryon ansatz reads 

\i;n) = \qnf\T.)\n)\x). (21) 

In the mean field approximation we demand the total energy functional E = {iph\ '■ H : \iph), 
where : if : is the normal ordered Hamiltonian of the models defined by (P, to be stationary 
with respect to variations of u{r), v{r), 0"(r), 0(r) and x(r). Of course, the meson wave 
functions may equivalently be determined by performing the variations with respect to the 
coherent state amplitudes C,{k), T]{k) and K{k). The variations with respect to the functions 
of r lead to a set of differential equations. For appropriate choices of the coupling constants, 
soliton solutions of those equations are obtained with three quarks absolutely confined (in 
CDM) gig or just bound (in LSM) [|,|Tg. 



The solitons described by the hedgehog state, \ijjh) cannot represent physical baryons be- 
cause they are not eigenstates of angular momentum or isospin. In addition, (PT| ) represents 
a localized object and therefore the translational symmetry of the model Hamiltonians is 
also broken in such states. In particular they contain spurious center-of-mass components 
which contribute to the energy and to other observables. 

States with good spin and isospin can be obtained from liph) by means of the Peierls- 
Yoccoz projection. The hedgehog only contains states with J = T and therefore, due to such 



space-isospace correlation, a single projection, either in spin or in isospin, is needed |[12| , |22 . 



We choose to project onto isospin. A baryon with isospin T, spin J = T and projection 
quantum numbers t and s (for isospin and spin, respectively) is given by 

|T,t;J = T, .) =(-lf+*pf_J^,), (22) 

where Vf_g is the isospin operator 

vL = ^^Jdnvl\n)R{n). (23) 

In this expression, R{Q) stands for the rotation operator in isospin space, V are the Wigner 
matrices and the integration is performed over all orientations Q (which represents the set 
of three Euler angles in isospin space). In the following we consider s = —t = M and use 
the shorthand notation Pjm = "Pm m- 

On the other hand, a model state representing a baryon at rest can be obtained by ap- 
plying a Peierls-Yoccoz projection onto linear momentum zero to the state already projected 
onto isospin (and angular momentum). The Peierls-Yoccoz linear momentum projector is 
given by 

P, = ^Jdae'^-^U{a), (24) 

where U{a) is the translation operator. A nucleon at rest is therefore represented by the 
model state 

J = T = ^,M,q = 0^ = P^=oPjM \A) = Pg=o \^jm) . (25) 

For q = 0, the isospin-angular momentum projector operator and the linear momentum 
projector operator commute, but this is no longer the case for q y^ 



In order to include recoil effects in the calculations, in principle one should boost [13 1 the 



zero momentum states (^51) , but the technical difficulties associated with boosting prevent, 
in practice, the use of such a procedure. However, at least for small linear momentum 
q, we may approximate the boost operation by the Peierls-Yoccoz projection onto linear 



momentum q [|TT],E3. Thus, the model state representing a physical baryon of angular 



momentum and isospin J and linear momentum q is 

\J = T,M,q)r^P,\i,jM) . (26) 

Proper normalization of the projected state requires the inclusion of kinematical normaliza- 
tion factors. For example a nucleon with four-momentum q, \N{q)), is described by 



\^N^{P^^j,,\P^ijj^,) 



(27) 



where J = \-i rriN is the nucleon mass and E = J q^ + m% its energy. 

Before presenting the formalism to compute the amplitudes for the electroproduction of 
the delta, one should briefly mention how the radial profiles M(r), f (r), cr(r), 0(r) and %('") 
are determined. They may be determined in the so-called "variation-before-projection" 



(VBP) method, and, in that case, the stationarity of the mean field energy is required. 
A better approach (even if much more demanding numerically) is the VAP ("variation- 
after-projection") method, where the energy functional to be minimized is the expectation 
value of the normal ordered Hamiltonian in the projected state \iP.]m)- In this procedure, 
which we followed, one obtains different field radial profiles UB^r), VB{r), (psir), crsir), 
XB^r) (and coherent state amplitudes ^^(fc), rjB{k), nsik)) for the nucleon {B -^ N) and 
for the delta [B -^ A). The results presented Section ^ use the VAP method for the 
angular momentum projection and the approximate VAP method for the linear momentum 
projection as described in [rT|. Unless otherwise stated, the coupling constants are ^f = 5 in 
the LSM |],|T2|,0 and G = 0.2 GeV in the CDM fgjgg], for which nucleon properties are 
well described. These values, together with the above mentioned values for meson masses 
and pion decay constant, will be referred to as the standard parameters. 



III. MULTIPOLE AMPLITUDES 



The N-A electromagnetic transverse helicity amplitudes PB| are defined by 



A 



(m) 



- / (A -A:A;a 



and the scalar helicity amplitude by 

Sx = —F- / (A - A; fcA 
V2v/2^\ 2 



J(0): 



J°(0) : 



A^-A 
2 



/x;fc 



N 



N -X-kN 



(28) 



(29) 



where J^ is the electromagnetic current density operator, e^, i = 0, ±1 are the photon 
polarization vectors (eo is chosen along the direction of the photon motion), and kw is 
the magnitude of the photon three-momentum at the photon point [Q. Because of gauge 
invariance, the other amplitude — longitudinal amplitude — is just the scalar amplitude 
multiplied by the kinematical factor ui/k. The values for A and /i are usually chosen as 
A = 1/2, 3/2, /i = 1 for the transverse amplitudes and A = 1/2 for the scalar amplitude. If 
the linear momentum projection of the model states is skipped (i.e. no recoil corrections). 



expressions (|28|) and (p9D reduce to those usually presented in the literature [see e.g. Eq. 
(7) of ref. 0] using the procedure described in [^ . 

Replacing the baryon states above by their model representations in Eq. (pT]), and noting 
that U^{r).r{Q)U{r) = r{r), we get (for A = 1/2, /i = 1) 



A 



1/2 



-N, 



S 



1/2 



V2 



NA 



N, 



f d^ad^re-''^'' /a— U^{x - l)a] : e-y ■ J{r) : U{xa 

U^[{x-l)a] : J%r) : U{xa) 



2 2 



d-'ad-'re""'' ( A 



11 



22 



A^- 



.11 
22 



(30) 
(31) 



where x is the fraction of the photon momentum carried by the delta. In this way, the 
parameter x identifies the reference frame used in the calculations: x = corresponds to 
the delta rest frame which is mostly used in the literature. The factor Nna contains all 
kinematical factors as well as the projected states normalization terms [see (p7[) ], and is 
given by 
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Nna = 



' E]\[ / E/\ 



(27r)3 v/2A^ V rriN V "^a ^Fi/2[{x - l)q]Fy2{xqy 



with 



Fi/2iq) 



N 



11 
22 



AT 



11 
22 



(32) 



(33) 



and similarly for the delta F factor. These factors only depend on the magnitude \q\. 

The multipole N — A transition amplitudes are usually extracted from the helicity am- 
plitudes above making a multipole expansion of the electromagnetic field. For nucleon and 
delta model states which are eigenstates of the angular momentum and parity operators, 
only the magnetic dipole and the electric and scalar quadrupole terms contribute to the 
transition (see Appendix for details). The linear momentum projection in our approach 
affects the rotational symmetry of the baryon states, and the relevant multipoles are not 
automatically selected. Instead, one has to explicitly remove the spurious terms in the 
multipole expansion of the exponential in equations (|3D|) and (^TJ), which implies restrict- 
ing the momentum transfer q to low values. In that case, the rotational symmetry of the 
model states is almost preserved even if the linear momentum projection is performed. The 
multipole amphtudes are then 



M^\(l) = -^NNAJd^ad^rniqr) 

U^[{x-l)a] : [f X J(r)]^ : U{xa) 



A 



1 1 



2 2 



N 



1 1 



M^'iq) 



VTo 



TT 



Nna f d^ad^r 

U^ix - l)a] : [V X J2(gr)r^2(^)] ■ J{r) : U{xa) 



2 2 



2 



(34) 



(35) 



M^'(g) 



VW^ Nna Jd^ad'^r j2 (gr ) F20 (r*) 

U^[{x-l)a] : J^{r) : U{xa) 



A 



1 1 



2 2 



N 



1 1 

2 2 



(36) 



where YJ-i are the vector spherical harmonics, ji{x) are the spherical Bessel functions and the 
index 1 in the Ml operator denotes component +1 in the spherical basis. It is worth noticing 
that formulas (|3^ ) -(|36| ) differ from those used when no recoil corrections are considered (see 
e.g. eqs. (10)-(12) of P) by the integration over a and by the presence of the translation 
operations. Had we inserted 6'^ (a) in (^)- (|3^) and integrated over a, the expressions for the 
multipole amplitudes when no recoil effects are considered (see, amongst others, [|6|, p!6| , p!8| , p!9| ) 
would be obtained. 

The electric quadrupole amplitude involves the operator 



0E2{q) = - fdM"^ X j2{qr)Yl2{r)] ■ J{r) 
q J 



(37) 



which, using the properties of the vector spherical harmonics and integration by parts, can 
be written as 



0E2(g) = ^ ^Jd'r^[rUqrW2iir)J'ir) 

-^ J d^rUqr) Y^^ir) r ■ J{r), (38) 

where we used the electric current conservation condition to simphfy the first term. The sec- 
ond term gives a neghgible correction to the E2 amphtude in the low momentum regime and 
can be dropped . Other technical aspects of the calculation of the multipole amplitudes 
are provided in the appendices ^ and 0. 

IV. RESULTS AND DISCUSSION 

The ratios E2/M1 and C2/M1 for the delta electroproduction are related to the multi- 
poles O-® through 

E2 _ 1 M^2 

mT~3M^ ^ ^ 

C2_ ^ J_M^ 

Ml 2V2M*^i' ^ ^ 

These ratios (EMR and CMR, respectively) are equal for \q\ ^ 0, a limit which is never 
met since, even at the photon point, a finite Iq^l is needed for the transition to take place. 

In most calculations reported in the literature the transition amplitudes are computed in 
the rest frame of the A. Such a choice corresponds to x = in the expressions of Section |T|. 
The nucleon four-momentum {E^, —q) and the photon four-momentum {u, q) completely 
specify the kinematics and the (invariant) photon virtuality, Q^ = — g^ is the appropriate 
quantity in terms of which the electroproduction amplitudes should be expressed. In the A 
reference frame, 

|qr=H + "^^ + ^'V-< (41) 



and 



UJ 



2m/ 



mi-m%- Q' 
2mA ' 



Figure |I| shows the results for the quadrupole electric to dipole magnetic ratio, in the 
LSM and CDM, for standard parameter sets in both models, in the rest frame of the delta. 
Figure ^ displays the quadrupole coulomb to dipole magnetic ratio as a function of —Q^. 
The first conclusion to be drawn is the compatibility of the model predictions with the 
data, namely the negative signs for both ratios. From the theoretical point of view we 
don't find any sign of the up and down behaviour of the data points. Another interesting 
conclusion is the small effect of the recoil corrections in EMR and CMR. Recoil corrections 



enhance the nucleon magnetic moments |Tl|,p^ and nucleon magnetic form factors [|14|. Such 



enhancement is also found in the nucleon-delta magnetic transition as it is shown in Figure 
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H, improving the comparison with experimental values, but the effect, in the present case, is 
smaller than for the nucleon. 

A similar enhancement turns out to show up in the quadrupole electric and coulomb 
multipoles, and altogether no sensible modification appears in EMR and CMR. In the CDM 
the modification of M^^ due to a better treatment of the kinematics of the nucleon and the 
delta is not enough to achieve a better comparison with the data. The comparison with the 
data of this observable favours the model with large number of pions in the cloud. The big 
slope of the theoretical CMR in the CDM is due to the small value predicted for M ^^^ in 
this model. 

The values of the ratios at the photon point (delta photoproduction) are —2.56% (LSM) 
and -2.54% (CDM) for the CMR and -2.11% (LSM), -1.85% (CDM) for the EMR. These val- 
ues are compatible (although slightly smaller, in the case of the EMR) with the experimental 
value —2.5 ± 0.5% estimated for EMR by the Particle Data Group p5| . 

It is not our purpose to find fittings of model parameters that better reproduce the 
experimental results (model parameters were fixed in the nucleon sector of the models). 
Nevertheless it is interesting to analyze the dependence of the results with model parameters 
namely the coupling constants. As stated before, the physical window for G in the CDM is 
relatively narrow and the resulting radial wave functions are very much similar throughout 
that physical range. The LSM, on the other hand, provides a larger range and a large 
variety of radial wave functions. The results are summarized in figures ^j-^ for three values 
of the coupling constant in the LSM: g = 4.5 (weak coupling, weak pion cloud), g = 5.0 
(intermediate coupling, standard parameter) and g = 5.5 (strong coupling, strong pion 
cloud). The graphs correspond to the calculation with recoil corrections. As Figure ^ 
reveals, the EMR remains impressively unchanged with —Q^. The CMR (Figure ^) is 
affected specially for large —Q^. The effect on Ml is shown in Figure ^ The multipole C2 
results from pion contribution alone, whereas Ml receives contributions from both pions 
and quarks. The stronger pion cloud enhances more C2 than Ml resulting in a larger (in 
absolute value) CMR for higher coupling constant. The same trends were also found in the 
CDM (but with even smaller variations with G). 

Finally we address the problem of the reference frame. In principle, the theoretical 
amplitudes should not be dependent on the particular choice of the reference frame. However, 
due to the lack of translational invariance of the model baryon states (even when recoil 
corrections are taken into account), that is not the case. Nevertheless, no dramatic changes 
in the results, as a consequence of the different choice of reference frame, are supposed to 
occur. In Figure ^ we present the Ml multipole amplitude and the CMR ratio for the 
LSM {g = 5.0) and for three values of the parameter x which is the fraction of the photon 
momentum carried by the A: x = 0, x = 0.5 and x = 0.7 (EMR follows the trend of CMR). 
The curve x = 0.5 would correspond to the Breit frame if nucleon and delta were degenerate. 
The major differences (indicating lack of covariance) come up at large values of —Q^ as one 
would anticipate. Indeed, unlike the correct description of baryon motion through Lorentz 
boosts of zero momentum eigenstates as mentioned in Section |T|, our approximate treatment 
is not relativistic and, therefore, more reliable for small and intermediate linear momenta. 
The region spanned by the curves in |^ gives an idea of the "theoretical uncertainty" of the 
model predictions. 
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V. CONCLUSIONS 

In this paper we addressed the question of the delta electroproduction amphtudes in 
the framework of two chiral effective models of the nucleon with meson and quark degrees 
of freedom. Although the predictions for the ratios E2/M1 and C2/M1 are compatible 
with data in both models, the amplitudes are better reproduced in the LSM thus favouring 
a picture of the nucleon and the delta with a stronger pion cloud. Recoil corrections of 
the baryons were taken into account in this study but no dramatic change was actually 
found with respect to the calculation with just angular momentum projection from the 
hedgehog. This is different from the modifications occurring in nucleon form factors where 
larger effects were found when the center of mass motion spurious components are removed 
from the baryon wavef unctions. Strong fluctuations on C2/M1, as seen in the experimental 
data, are not observed in the present approach. 

Acknowledgements: We thank B. GoUi, M. Rosina, S. Sirca, A. Silva and D. Urbano for 
useful discussions. Financial support from Fundagao para a Ciencia e Tecnologia, Portu- 
gal, project PRAXIS/PCEX/P/FIS/6/96 is acknowledged. One of the authors (LA) also 
acknowledges the financial support to his PhD program from PRODEP Project No. 185/007. 

APPENDIX A: DEFINITION OF THE MAGNETIC AND ELECTRIC 

MULTIPOLE AMPLITUDES 

The operator involved in the calculation of the transverse helicity amplitudes is 

d{r) = e-^-^'ei ■ J{r). (Al) 

Choosing the z-axis in the direction of q, the expansion of the exponential reads 

e--^- = v^5:(-z)'v/2lTTj/(gr)r,o(r), (A2) 

where ji{x) are the spherical Bessel functions and Yim{r) are the spherical harmonics. Now, 
the product Yi^{r)ei, can be cast in terms of the vector spherical harmonics as 

Yi^{r)e, = J2{lfi; lu\jm) Y]}ir), (A3) 

jm 

where {lfi;l'fi'\jm) are Clebsch-Gordan coefficients. Equation ( |AT| ) then reads 



0{r) = v^^(-0'v/27njKgr)(/0; ll|jl)l^],(f) ■ Jir). (A4) 

ij 

Since only the terms with j = I, j = I ±1 contribute, we may write 



d{r) = v^J: v/2m(-0^ Ln{qr)YUr) ■ J{r) 



+ 1 



^ ^ ^-jL-i{qr)Yl,^,{r) • J{r) - \U—JL+i{qr)Y LL+i{r) ■ J{r) 



2L + r '^ ' ^"^"'^ ' ^ ' ^ 2L + r 

1, 



^5: v/2m(-0^ \ -jL{qr)YUr) + -Vx [jL{qr)YUr) \ ■ J{r) (A5) 



L 
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(see [^ for more details). The two terms inside the curly braces in this expression are 
respectively the electromagnetic field {L, l)-magnetic and (L, l)-electric multipoles 



A[%ir) 



jL{qr)Yf^{r 



-k^"" 



jL{qr)Y^Ur 



(A6) 
(A7) 



As it is shown for instance in |3^, the scalar products of the L-th order field multipoles with 
any vector (like the current density operator) form the irreducible components of rank-L 
operators, with parity (—1)^"'"^ and (—1)^ for the magnetic and electric multipoles, respec- 
tively. In a transition between states with angular momentum Jj = 1/2 and J/ = 3/2 and 
positive parity, only (L = l)-magnetic and (L = 2)-electric multipoles may contribute, so 
that the operator O may be replaced by 



0'(r) 



/67iA\T'ir) - VlO^Ail\r] 



dM)i 



J{r). 



(A8) 



The Ml and E2 amplitudes are, respectively, the matrix elements of the first and second 
terms on the right hand side of this equation. To make the correspondence with equation (|34D 
we note that 



3z 



iVQ^AYC\r) ■ J{r) = ^Mqr) ^(1/i; lu\ll)f^Mr) 



V2 
3 



flU 



ji(gr) [r x J{r)]^ 



(A9) 



where use was made of the definition of the vector spherical harmonics and of the expression 
of the spherical components of the vector product of two vectors. The scalar amplitude can 
be derived in a similar fashion. 



APPENDIX B: TRANSITION OVERLAP 

The transition overlap of two (not necessarily the same) hedgehog baryons, defined as 

N{a,Q) = {^'^\U{a)RmiJh), (Bl) 

is a recurring function in calculations involving isospin and linear momentum projected 



states |[TT|j2J] . It is the following product of quark and meson overlaps: 

N{a, n) = N^{a, fi) N^{a, n)N^{a, fi) N^{a, fi). 
An explicit form for this function can be derived following the calculation of the norm 



overlaps in ref . [^ . The isoscalar meson overlaps do not depend on Q or on the orientation 
of a, and taking advantage of the properties of the coherent states (for the sigma field for 
instance), one obtains 



N^{a, n) = n^{a) = exp <^ g^{a) - n dk k^ uj^{k) [r]' (k) + r]^{k) 
where we introduced the functions 



(B2) 
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g^{a) = 27r / dk k^ uj^{k)ji{ka)ri'{k)ri{k) . 



(B3) 



A similar function should also be defined for the chromodielectric field. The quark overlap 
is also readily computed, because the spatial part of the quark wavefunctions is invariant 
under isospin rotations and the spin-isospin part is invariant under space translations. It is 
given by 



2 r°° 

dk k^joika) [u'{k)u{k) + v\k)v{k)] 



nq[a) = - 

n JO 

\r fn\ ^ " + ^ 

J\q\}l) = COS — COS . 

In these equations, u and v are Fourier transforms of the quark profiles, given by 

/•oo 

u(k) = / dr r"^ jo(kr)u{r) 



dr r"^ ii{kr)v{r). 



v{k) = 
For the pion field overlap, we get 

iV.(a,f]) =exp{-7r|^"rffcfcV(A:) [^^k] + eik)]'^ exp |1 [^J(a) + (^^"(a)] Tri?(f])} 



(B4) 
(B5) 

(B6) 

(B7) 
(B8) 



exp{g^{a) diRij{Q)dj} , 
with gf defined as for the sigma (see Eq. p3|) : 

gUa) = 27r / dk ^ uMji{ka)e{knk). 



(B9) 



(BIO) 



In the reference frame with the z-axis along the axis of the rotation R, the quantity 
diRij{Q)dj does not depend on the azimuthal angle of vector a and that could be exploited 



in order to simplify some integrations [11|. However, the orientation of the z-axis has already 
been fixed along the direction of the photon momentum. Therefore, in all integrations over 
a the transformation a -^ T^^a is made, where T is the rotation that aligns the z-axis with 
the axis of the rotation i?(f2), and again advantage can be taken from the above mentioned 
independence of the azimuthal angle. One gets 



N^{T ^a, n) = n^{a)J\f^{a, s, Vl) , 
where s = cos 6a is the cosine of the polar angle of a, and 

n^{a) = exp (^q («) - t^ T dk k^uj^{k) k'^{k) + ^\k) 



AC(a, s, fi) = exp < 2z{a, s) 



2 



2 H 2<^ + l 
COS — COS 



P2{s) being the Legendre polynomial of second degree. 



(Bll) 

(B12) 
(B13) 
(B14) 
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APPENDIX C: CALCULATION OF MATRIX ELEMENTS WITH PROJECTED 

STATES 

Here we present some details regarding the calculation of matrix elements of opera- 
tors between isospin and linear momentum projected states, using the C2 amplitude - see 
Eq. (|36| ) - as an example. The other amplitudes can be obtained in a similar way (details 
ini3). 

The electromagnetic current density for the two effective theories considered in this 
paper, derived using Noether's theorem, reads 

3 



J'ir) = j:Ac)ir)Y (q + 2^f^) ^(c)W + [^(r) x d^n{r)],, (CI) 



where c is a quark index and the cross product in the second term is in isospin space. The 
current density is the sum of an isoscalar operator 

S'{r) = lJ:^(.){r)Y^^,){r), (C2) 

which, because of isospin conservation, can not contribute to A^-A matrix elements, and the 
zeroth component of an isovector operator 

y^A^) = ^E^(c)W7^r(^V(c)(r) + [7f(r) x d^^[r\. (C3) 

^ c=l 

The components of isovector operators commute in a well defined manner with the isospin- 



space rotations involved in the isospin projectors (see |3^, for instance) and one can show 
that, regarding the expression of the C2 amplitude, 

i 

with 



"' = vl(i' + !''•-* 



f4>- («) 



We can then write the C2 amplitude as 



M''\q) = -^N^AT.^tJd'a J dQ P%+,) _y,{n) F,(a,fi), (C6) 

with 

Ft{a,n) = [ £r Uqr)Y2oir){M^)\U^[ix - ^)a]ViAr)Rin)U{xa)\MN)). (C7) 



In this expression, \iph{N)) and \iph{^)) represent the nucleon and the delta hedgehogs. 

The quark component of the isovector part of the charge density in Eq. (^^) cannot 
contribute here since the C2 is a matrix element of an L = 2 operator between s-wave 
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quark states. We are then left only with the pion contribution to the charge density, and 
we expand the pion field 7f(r) and its canonical conjugate Pni^) in plane waves: 

1 /■ .,- 1 



TTi{r) 



PUr) 



d'k- 



(2-)3/^ J ^^^^) 



ai{k)e"^"^ + al{k)e 



-ikr 



f27r)3/2 



3,../^^(^) r 



d-'k 



kr f/ 



ai{k)e"^"^ - al{k)e 



^—ikr 



(C8) 



(C9) 



The translations and/or isorotations of coherent states still yield coherent states (the trans- 
formed states are still eigenstates of the annihilation operators) with shifted and/or isoro- 
tated amplitudes. Indeed, for the pion field one has 



ai{k)U{xa)R{n)\i,f,{N)) 



luJn{k) 



-ixk-a 



.(N) 



R,,{n) e/ >{k) U{xa)R{n)\MN)), (CIO) 



(^.(A)| U^[ix-l)a] a,(fc) = .j^M e^i^-^)^''^ ^f\k) {M^)\ U^[{x-l)a]. (Cll) 

Taking now advantage of the hedgehog shape of the pion coherent state amplitude, 
Eq. (|18D, we can write the cartesian components of the functions -Fj(a, VL) in Eq. ( |C7|) as 

F,(a,fi) = -^iV(a,fi)y'rfV j2(gr)y2o(r')e.,fci?fcz(fi)(r-)Kr+),«(a,r,a-f), (C12) 



where 



r+ = r + (1 — x)a, 
r_ = r — xa, 

., A^(r_)5A(r+) + AA(r+)5jv(r. 



r_r I 



AJr) 



dkk^j,{kr)Uk), 
Bnir) = I dk k^u^ik)j,ikr)Uk) 



(C13) 
(CM) 

(C15) 

(C16) 
(C17) 



and 



N{a,n) = {MA)\U{a)R{n)\MN)) (C18) 

is the transition overlap of the A^ and the A intrinsic hedgehogs, derived in Appendix |B|. 
Expanding the products (r_)i (r+)j, the function ( |C12[ ) unfolds in three terms, 

F,(a, n) = Fl°\a, fi) + F}^\a, fi) + F^ia, fi) , (C19) 

with 

Ff)(a,n) = ^N{a,n) a^ x{l - x) ei,kRki{n)djdi J d^r j2{qr)Y2o{r)a{a,r,a ■ r) (C20) 

F^'\a,n) = -^N{a,n)ae^,kRki{^) 

ZTT 

X d^rrJ2{qr)Y2o{r)[{l — x)djfi — xdifj]a{a,r,a-r) (C21) 

FPia,n) = -^Nia,n)e,jkRkiin) fd'rr^j2iqr)Y2oir)rifja{a,r,a-r). (C22) 

ZTT J 
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Let us focus on F-- . In doing the integration over r we are not allowed to pick any 
particularly convenient orientation for the vector a, whose components are also integration 
variables. Still, it is possible to perform the integration over the azimuthal angle of r 
analytically, yielding 

FfV,^) = 1^^ x{l-x) a'So2ia) e,,kRkim N{a,n) (3a^ - l) afa^ (C23) 

where the following functions were introduced 

SUa)= drr^+^Uqr) du Piiu)a{a,r,u), (C24) 

Jo J-i 

Pi{u) are the Legendre polynomials and u = cos 6r is the cosine of the polar angle of r. 
Further analytical refinations of these F functions are not possible because the function 
a{a, r, u) depends on the radial profiles of the meson fields, which are known only numerically. 
In order to compute the contribution of the function F^^^ to M'^'^, we replace its spherical 
components in Eq. (|C6|) . Making the transformation a -^ T~^a (see Appendix ^) the 
dependence of the integrand on the azimuthal orientation of the vector a is restricted to 
the terms ajdi and d'^djdi, allowing us to perform the integration over (pa analytically. The 
contribution of F^^^ to the Coulomb amplitude can then be written as 



^C2(0)/_N V^IO^ 



(q) = -j^ N^A E c* /" da a' f_^ ds J rffiP%+,) _y,{n) 4'\a, s, Q), (C25) 



(27v 
where s = cos 6a is the cosine of the polar angle of a and 



r2n 



/f)= f^d<j)aF,^'\T-'a,n). (C26) 

Jo 

The integration over the azimuthal angle of a can now be performed analytically, as the 
dependence on this variable only appears in the vector component products such as djdi, 
etc. The cartesian components of the isovector function I^^^ are then 

Jj- (a, s, fi) = i\ —r x{l — x) 0^5*02(0.) n{a)N{a, s, fi) 

3P^{s) T,{n) - ^ [Po{s) - P2{s)] x,{n) - P2{s)w,{n) 

+^ [7Po{s) - lOP2{s) + 3Pi{s)] [2Y,{n) + X^m 
3 



where 



+ ^ [P2{s) - P^is)] [2Zi{Q) + 2Vi{Q) + Wi{Q) + t/,(fi)]| , (C27) 



N{a,s,Q) = N{T-^a,Q) (C28) 



(see the Appendix |B|) and the following isovectors, functions of the Euler angles alone, have 
been introduced 

17 



Ti{n) = e,,kRki{n)T,j{n)Tsi{n)Ti^{n) 
Viin) = e^,kRk3{^)Tsjin)T33in) 

Xi{U) = eijkRkji^) 

z,in) = e,,kRu{n)niin)T,,in). 

In these definitions, T,-, are the matrix elements of the transformation T. We now note 



-«j 



that RkiT^i = Tf,^ Tsk'Rk'iTi^ = T^k, because Tgk'Rk'iTi 
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^s3 



5s3 is the rotation 



matrix transformed to the frame in which the z-axis is along the axis of rotation. Then, 
we immediately obtain Ti = 0, Wi = and Zi = eiskT^kTss. Finally, we reconstruct the 
spherical components of J^^-* , replace the result in Eq. ( |C25| ) and, after the contraction with 
the Wigner matrices V, the integration over the Euler angles can proceed analytically, along 
the lines followed in ref. ITTIl. The final result is 



M'^^^'"' (g) 



IT 



6\/2 



x{l - x) N, 



NA 



da a^ n{a) So2{a) [Tio(a) - T^ia) - Tso^a) + T32(a)] 



(C29) 



with Tki defined by 



niia) = r ds —I,{z)Pi{s) 
J-i z 



(C30) 



where Ik are the modified Bessel functions of the first kind, depending on the parameter z 
introduced in ( |B14| ). Using a similar procedure we computed the two remaining contributions 
to the scalar amplitude. The whole method can be applied to obtain formulas for the 
magnetic dipole amplitude fS^ . 
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FIG. 1. Quadrupole electric to dipole magnetic ratio in the LSM and the CDM, as a function 
of —Q^. The sohd (dashed) hnes show the results with (without) recoil effects for the standard 



parameter set and in the rest frame of the A. Experimental data was taken from |27,2^ (o) 
and mM (•)• 
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FIG. 2. Quadrupole Coulomb to dipole magnetic ratio. Conventions, parameters and reference 
frame as in Fig. ||. Experimental data taken from |27,3^ (o) and [30,3^ (•). 
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FIG. 3. Dipole magnetic amplitude Ml (in units of 10 ^ GeV ^'^) as a function of —Q^, for 
the LSM and the CDM. Conventions, parameters and reference frame as in Fig. |l[ Experimental 



data taken from |34| 
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FIG. 4. EMR in the LSM for three values of the coupling constant. 
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FIG. 5. CMR in the LSM for three values of the couphng constant. 
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FIG. 6. Ml in the LSM for three values of the coupling constant. 
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FIG. 7. Ml (left) and CMR (right) in the LSM, computed in different inertial frames using the 
standard parameter set in the LSM. The parameter x denotes the fraction of the photon momentum 
carried by the A. 
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